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ABSTRACT 

Aims.  We  investigate  the  structure  and  shape  of  the  photospheric  and  molecular  layers  of  the  atmospheres  of  four  Mira  variables. 
Methods.  We  obtained  near-infrared  K-band  spectro-interferometric  observations  of  the  Mira  variables  R  Cnc,  X  Hya,  W  Vel,  and 
RW  Vel  with  a  spectral  resolution  of  about  1500  using  the  AMBER  instrument  at  the  VLTI.  We  obtained  concurrent  JHKL  photometry 
using  the  the  Mk  II  instrument  at  the  SAAO. 

Results.  The  Mira  stars  in  our  sample  are  found  to  have  wavelength-dependent  visibility  values  that  are  consistent  with  earlier 
low-resolution  AMBER  observations  of  S  Ori  and  with  the  predictions  of  dynamic  model  atmosphere  series  based  on  self-excited 
pulsation  models.  The  corresponding  wavelength-dependent  uniform  disk  (UD)  diameters  show  a  minimum  near  the  near-continuum 
bandpass  at  2.25  //m.  They  then  increase  by  up  to  30%  toward  the  H2O  band  at  2.0 /in 1  and  by  up  to  70%  at  the  CO  bandheads  between 
2.29  fj m  and  2.48  /r m.  The  dynamic  model  atmosphere  series  show  a  consistent  wavelength-dependence,  and  their  parameters  such 
as  the  visual  phase,  effective  temperature,  and  distances  are  consistent  with  independent  estimates.  The  closure  phases  have  signif¬ 
icantly  wavelength-dependent  and  non-zero  values  at  all  wavelengths  indicating  deviations  from  point  symmetry.  For  example,  the 
R  Cnc  closure  phase  is  1 10°  ±  4°  in  the  2.0  //m  H2O  band,  corresponding  for  instance  to  an  additional  unresolved  spot  contributing  3% 
of  the  total  flux  at  a  separation  of  ~4  mas. 

Conclusions.  Our  observations  are  consistent  with  the  predictions  of  the  latest  dynamic  model  atmosphere  series  based  on  self-excited 
pulsation  models.  The  wavelength-dependent  radius  variations  are  interpreted  as  the  effect  of  molecular  layers  lying  above  the  pho¬ 
tosphere.  The  wavelength-dependent  closure  phase  values  are  indicative  of  deviations  from  point  symmetry  at  all  wavelengths,  thus 
a  complex  non-spherical  stratification  of  the  extended  atmosphere.  In  particular,  the  significant  deviation  from  point  symmetry  in  the 
H20  band  is  interpreted  as  a  signature  on  large  scales  (there  being  a  few  across  the  stellar  disk)  of  inhomogeneities  or  clumps  in 
the  water  vapor  layer.  The  observed  inhomogeneities  might  possibly  be  caused  by  pulsation-  and  shock-induced  chaotic  motion  in 
the  extended  atmosphere. 

Key  words,  techniques:  interferometric  -  techniques:  photometric  -  stars:  AGB  and  post-AGB  -  stars:  atmospheres  - 
stars:  fundamental  parameters  -  stars:  mass-loss 


1.  Introduction 

Mira  stars  are  long-period,  large-amplitude  variable  stars  on  the 
asymptotic  giant  branch  (AGB).  Mass  loss  becomes  increasingly 
important  toward  the  tip  of  the  AGB,  before  the  star  evolves 
to  the  planetary  nebula  (PN)  phase,  where  a  great  diversity  of 
morphologies  is  seen.  Contemporary  astrophysical  questions  re¬ 
lated  to  the  study  of  Mira  variables  include  the  pursuit  of  the 
detailed  mass-loss  mechanism  on  the  AGB,  including  the  effects 
of  pulsation  and  shock  fronts  on  the  structure  and  morphology 


*  Based  on  observations  made  with  the  VLT  Interferometer  (VLTI) 
at  Paranal  Observatory  under  program  ID  082.D-0723. 

**  Figure  2  is  available  in  electronic  form  at 
http : //www . aanda . or g 


of  the  extended  atmosphere,  and  of  the  mechanism  shaping  the 
observed  PN  morphologies.  Near-infrared  interferometry  has 
traditionally  been  used  to  characterize  AGB  star  atmospheres. 
In  particular,  observations  using  the  IOTA  interferometer  have 
uncovered  the  wavelength-dependence  of  Mira  star  diameters 
using  a  few  bandpasses  with  spectral  resolutions  of  up  to  A/ AT  ~ 
25,  providing  observational  evidence  of  molecular  layers  lying 
outside  the  photospheric  layers  (e.g.;  Mennesson  et  al.  2002; 
Perrin  et  al.  2004).  These  molecular  layers  have  also  been 
present  in  theoretical  dynamic  model  atmospheres  (Hofmann 
et  al.  1998;  Ireland  et  al.  2004a,b).  H- band  interferometry  at  the 
IOTA  interferometer  in  the  broad  band  (Ragland  et  al.  2006), 
as  well  as  in  three  filters  with  4 /AT  ~  15  (Ragland  et  al. 
2008;  Pluzhnik  et  al.  2009)  have  revealed  non-zero  closure 
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Table  1.  Log  of  our  VLTI/ AMBER  observations. 


Target 

Date 

JD 

Mode 

Baseline 

m 

PA 

deg 

R  Cnc 

2008-12-29 

4830 

MR23 

15.9/31.8/47.7 

-74 

2008-12-30 

4831 

MR21 

16.0/32.0/48.0 

-73 

X  Hya 

2008-12-29 

4830 

MR23 

16.0/31.9/47.9 

-71 

2008-12-30 

4831 

MR21 

15.5/30.9/46.4 

-76 

W  Vel 

2008-12-29 

4830 

MR23 

15.3/30.6/45.9 

-84 

RW  Vel 

2009-03-01 

4892 

MR23 

15.1/30.1/45.2 

-86 

Notes.  The  Julian  date  (JD)  is  given  as  JD-2  450  000.  The  AMBER  in¬ 
strument  mode  is  MR23  (medium  resolution  mode  with  central 
wavelength  2.3  /mi  and  range  2. 12-2.47  pm)  or  MR21  (2.1 /.on, 
1.92-2. 26 pm).  The  baseline  is  the  projected  baseline  length  for  the 
AT  VLTI  baselines  used  (E0-G0/G0-H0/E0-H0)  with  projected  base¬ 
line  angle  PA. 


Table  2.  SAAO  JHKL  photometry. 


Target 

J 

mag 

H 

mag 

K 

mag 

L 

mag 

Ay 

mag 

mb  01 

mag 

R  Cnc 

0.866 

-0.155 

-0.652 

-1.063 

0.07 

2.55 

X  Hya 

2.375 

1.476 

0.982 

0.397 

0.18 

4.12 

W  Vel 

2.002 

0.987 

0.503 

-0.097 

0.42 

3.64 

Notes.  The  observations  were  obtained  on  2008  December  30 
(JD  2  454  831),  concurrently  with  the  AMBER  observations  in  Table  1 . 

phases  for  several  Mira  variables,  which  reflect  asymmetric 
brightness  distributions  of  the  photosphere  and/or  the  enve¬ 
lope  around  the  star.  Wittkowski  et  al.  (2008)  presented  the 
first  VLTI/ AMBER  near-infrared  spectro-interferometric  obser¬ 
vation  of  an  AGB  star,  providing  continuous  wavelength  cov¬ 
erage  from  1 .29 /./m  to  2.32 pm  with  a  spectral  resolution  of 
/1/A/l  ~  35.  The  data  showed  visibility  and  diameter  variations 
as  a  function  of  wavelength  that  generally  confirmed  the  predic¬ 
tions  of  dynamic  model  atmospheres,  where  the  diameter  varia¬ 
tions  can  be  understood  as  the  effects  from  atmospheric  molec¬ 
ular  layers  (most  importantly  FLO,  CO,  and  TiO).  Here,  we 
present  the  first  VLTI/ AMBER  observations  of  Mira  variables 
using  its  K  medium  resolution  modes  with  a  spectral  resolution 
of  4/ AT  ~  1500,  and  a  comparison  to  the  newly  available  CODEX 
dynamic  model-atmosphere  series  by  Ireland  et  al.  (2008,  2011), 
which  are  based  on  self-excited  pulsation  models. 

2.  Observations 

We  obtained  K- band  spectro-interferometry  with  a  resolution  of 
T/AT  ~  1500  of  the  Mira  variables  R  Cnc,  X  Hya,  W  Vel, 
and  RW  Vel  with  the  AMBER  instrument  (Petrov  et  al.  2007) 
of  the  VLTI.  Table  1  shows  the  details  of  our  observations.  We 
obtained  concurrent  near- infrared  JHKL  photometry  of  R  Cnc, 
X  Hya,  and  W  Vel  at  the  SAAO  Mk  II  instrument.  Table  2  lists 
the  obtained  photometry  and  the  estimated  bolometric  magni¬ 
tudes  using  the  procedures  outlined  in  Whitelock  et  al.  (2008) 
and  their  estimates  of  Ay  -  For  R  Cnc  and  X  Hya,  we  obtained  ad¬ 
ditional  AMBER  and  SAAO  observations  at  later  epochs/phases, 
as  well  as  coordinated  mid-infrared  VLTI/MIDI  observations 
and  VLBA  observations  of  the  SiO  and  H2O  maser  emission. 
The  complete  data  set  for  these  sources  will  be  discussed  in  a 
dedicated  forthcoming  paper. 

Raw  visibility  and  closure  phase  values  were  obtained  from 
the  AMBER  data  using  version  3.0  of  the  amdlib  data  reduc¬ 
tion  package  (Tatulli  et  al.  2007;  Chelli  et  al.  2009).  An  absolute 


wavelength  calibration  was  performed  by  correlating  the 
AMBER  flux  spectra  with  a  reference  spectrum  that  includes 
the  AMBER  transmission  curve,  the  telluric  spectrum  estimated 
with  ATRAN  (Lord  1992),  and  the  expected  stellar  spectrum  us¬ 
ing  the  BS  4432  spectrum  from  Langon  &  Wood  (2000),  which 
has  a  similar  spectral  type  as  our  calibrators  (K3-4  III).  Similarly, 
a  relative  flux  calibration  of  the  target  stars  was  performed  using 
the  calibrator  stars  and  the  BS  4432  spectrum.  Calibrated  visi¬ 
bility  spectra  were  obtained  by  using  an  average  of  two  trans¬ 
fer  function  measurements  taken  before  and  after  each  science 
target  observation.  The  calibrators  were  l  Hya  (for  R  Cnc  and 
X  Hya),  81  Gem  (for  R  Cnc),  31  Leo  (for  X  Hya),  GZ  Vel  (for 
W  Vel  and  RW  Vel),  and  y  Pyx  (for  W  Vel),  and  were  selected 
from  Borde  et  al.  (2002)  using  the  ESO  calibrator  tool  CalVin. 

3.  Results  and  comparison  to  model  atmospheres 

Figures  1  and  2  show  the  resulting  AMBER  flux,  visibility, 
UD  diameter,  and  closure  phase  data.  The  UD  diameter  curve  is 
obtained  by  fitting  a  UD  model  to  the  visibility  data  separately 
for  each  spectral  channel.  We  note  that  the  intensity  profiles  are 
expected  to  be  more  complex  than  a  simple  UD  model,  and  this 
approach  is  used  merely  to  readily  obtain  a  preliminary  estimate 
of  the  apparent  diameter.  Signatures  of  a  more  complex  model  as 
well  as  asymmetries  are  most  pronounced  near  and  beyond  the 
first  visibility  minimum,  hence  we  used  the  two  shorter  baselines 
only  for  UD  fits. 

The  AMBER  data  are  compared  to  the  P/M  dynamic  model 
atmosphere  series  (Ireland  et  al.  2004a, b),  as  well  as  the  newly 
available  CODEX  series  (Ireland  et  al.  2008,  2011).  Both  series 
are  based  on  self-excited  pulsation  models.  The  only  fit  parame¬ 
ter  for  each  model  atmosphere  is  the  Rosseland  angular  diameter. 
Figure  1  shows  the  result  of  the  best-fit  model  of  each  of  the  P/M 
and  CODEX  series.  Compared  to  the  P/M  series,  the  CODEX  se¬ 
ries  uses  the  opacity  sampling  method,  and  the  model  series  are 
available  for  a  wider  range  of  basic  parameters  of  the  underly¬ 
ing  hypothetical  non-pulsating  parent  star.  The  CODEX  intensity 
profile  was  tabulated  in  steps  of  0.0005  pm  (T/AT  =  40  000  at 
A  -  2  pm),  and  the  synthetic  visibility  values  averaged  over  the 
AMBER  spectral  channels  assuming  a  Gaussian-shaped  band¬ 
pass.  The  CODEX  series  include  the  series  o54  (parent  star  param¬ 
eters:  M  =  1.1  MQ,  L  =  5400 Lq,  R  =  216R0,  P  =  330  days), 
R52  (1.10Mo,  5200  L0,  209 R0,  307  d),  C5®  (1.35 M0,  5050 L0, 
29 1  Rb,  430  d),  and  C8 1  ( 1 .35  M0,  8 1 60  L0, 278  R0, 427  d),  com¬ 
pared  to  the  previous  series  P  (1.0Mo,  3470  L0,  241  RQ,  332  d) 
and  M  (1.2  Me,  3470  L0,  260 Re,  260  d).  Table  3  provides  an 
overview  of  the  measured  and  derived  quantities  of  the  four  tar¬ 
gets. 

The  targets  of  our  sample  have  a  visibility  function  with  a 
characteristic  shape1.  For  each  baseline,  it  exhibits  a  maximum 
near  2.25  pm  and  decreases  towards  both  shorter  and  longer 
wavelengths.  At  shorter  wavelengths,  it  shows  a  minimum  near 
2.0  pm.  At  longer  wavelengths,  between  2.3  pm  and  2.5  pm,  the 
visibility  function  displays  sharp  drops  at  the  locations  of  the 
CO  bandheads.  This  characteristic  shape  of  the  visibility  func¬ 
tion  of  the  four  targets  studied  here  is  consistent  with  low  res¬ 
olution  (T/AA  ~  35)  AMBER  data  of  the  Mira  variable  S  Ori 
(Wittkowski  et  al.  2008),  but  shows  additional  details  thanks  to 
the  increased  spectral  resolution.  Figure  1  displays  the  position 
of  the  FLO  (1.9-2.  lpm)  and  CO  bands  (2. 3-2. 5 pm,  but  where 
H2O  may  also  contribute). 

1  The  R  Cnc  visibility  reaches  into  the  second  lobe,  where  the  shape  is 
inverted. 
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Wavelength  (nm) 


Fig.  1.  R  Cnc  flux,  squared  visibility  amplitude,  closure  phase,  and 
UD  diameter  (from  top  to  bottom )  as  a  function  of  wavelength.  Data  of 
the  remaining  three  targets  are  shown  in  Fig.  2  of  the  electronic  version. 
The  wavelength  range  is  1.92-2.47  pm  (MR21  mode:  1.92-2.26  pm, 
MR23  mode:  range  2. 12-2.47  gm).  The  modes  overlap  in  the  region 
2. 12-2.26  pm,  where  the  two  lines  appear  as  one  thicker  line.  The  thick 
vertical  lines  denote  the  errors  averaged  over  three  wavelength  intervals. 
The  blue  lines  show  a  comparison  to  the  best  fitting  model  of  the  P/M 
model  atmosphere  series  (Ireland  et  al.  2004a, b),  and  the  red  lines  the 
best  fitting  model  of  the  new  CODEX  series  (Ireland  et  al.  2008,  2011). 
The  dashed  line  shows  a  uniform  disk  curve  fitted  to  the  2.25  pm  near¬ 
continuum  bandpass. 


The  corresponding  wavelength-dependent  uniform  disk 
(UD)  diameters  show  a  minimum  near  the  near-continuum  band¬ 
pass  at  2.25 pm.  They  increase  toward  2.0pm ,  where  the  FUO 
opacity  is  large,  by  up  to  ~30%.  The  CO  bandheads  between 
2.29 pm  and  2.48 pm  are  clearly  visible  for  all  targets  with  an 
increase  in  the  UD  diameter  by  up  to  ~70%. 


The  best-fit  models  of  the  P/M  and  CODEX  model  series  pre¬ 
dict  visibility  and  UD  curves  that  are  consistent  with  the  obser¬ 
vations.  Here,  the  newly  available  models  of  the  CODEX  series 
provide  a  closer  agreement  with  the  data  than  the  P/M  series,  at 
both  the  locations  of  the  CO  bandheads,  which  may  be  expected 
because  of  the  newly  introduced  opacity  sampling  method  and 
their  higher  spectral  resolution,  and  in  terms  of  the  overall  shape 
of  the  curve.  The  latter  may  be  explained  by  the  availability  of 
more  model  series  with  additional  parent  star  parameters  and 
with  an  increased  phase  coverage  per  series.  The  AMBER  flux 
curves  are  consistent  with  those  of  the  dynamic  model  atmo¬ 
sphere  series.  In  addition,  the  pulsation  phases  of  the  best-fit 
CODEX  models  are  consistent  with  the  observed  visual  phases, 
and  their  effective  temperature  is  consistent  with  that  derived 
from  the  fitted  angular  Rosseland-mean  diameter  and  the  bolo- 
metric  flux.  Furthermore,  the  distances  obtained  from  the  fitted 
angular  diameter  and  the  model  radius  are  consistent  with  those 
derived  from  the  measured  apparent  bolometric  magnitude  and 
the  model  luminosity,  as  well  as  with  the  period-luminosity  dis¬ 
tance  from  Whitelock  et  al.  (2008,  cf.  Table  3)  within  l-3<x. 

The  closure  phase  functions  of  our  targets  exhibit  significant 
wavelength-dependent  non-zero  values  at  all  wavelengths.  Non¬ 
zero  values  of  the  closure  phase  are  indicative  of  deviations  from 
point  symmetry,  here  along  the  projected  intensity  profile  onto 
the  orientation  of  the  position  angle  given  in  Table  1 .  The  consis¬ 
tency  of  the  visibility  amplitudes  with  spherical  models  and  the 
closure  phase  deviation  being  more  significant  for  targets  that 
are  well-resolved  indicates  that  the  deviation  from  point  sym¬ 
metry  originates  from  sub-structure  at  a  relatively  low  flux  level 
of  an  overall  spherical  intensity  distribution.  The  strongest  clo¬ 
sure  phase  signal  is  obtained  for  the  most  clearly  resolved  target 
R  Cnc.  It  differs  strongly  from  point- symmetry  in  the  HoO  band 
at  2.0pm  with  a  value  of  1 10°  +  4°.  At  the  continuum  bandpass 
at  2.25  pm.  the  closure  phase  value  is  closer  to  0  with  a  value 
of  45° +  5°.  Towards  longer  wavelengths  (CO  band),  the  aver¬ 
age  value  is  153°  ±  4°,  with  values  up  to  165°  in  the  bandheads, 
and  thus  close  to  a  180°  symmetric  value.  As  an  illustration,  the 
R  Cnc  closure  phase  value  in  the  HiO  band  can  be  caused  by 
the  addition  of  an  unresolved  (up  to  ~3  mas  diameter)  spot  to  a 
circular  disk,  contributing  up  to  ~3%  of  the  total  flux,  i.e.  cor¬ 
responding  to  the  level  of  the  visibility  departures  between  data 
and  spherical  models,  at  a  separation  of  ~4  mas. 

4.  Discussion 

The  four  Mira  variables  of  our  sample  exhibit  consistent  char¬ 
acteristic  wavelength  dependences  of  the  visibility  and  conse¬ 
quently  the  corresponding  uniform  disk  diameter  that  are  con¬ 
sistent  with  those  of  earlier  low-resolution  AMBER  data  of  the 
Mira  variable  S  Ori  and  the  predictions  of  the  P/M  and  CODEX 
dynamic  model  atmosphere  series.  Here,  the  newly  available 
CODEX  series  provides  a  closer  agreement  with  the  data  than 
the  earlier  P/M  series.  This  result  confirms  that  the  wavelength- 
dependent  angular  diameter  is  caused  by  the  atmospheric  molec¬ 
ular  layers,  here  most  importantly  H2O  and  CO,  as  they  are 
naturally  included  in  the  dynamic  model  atmosphere  series. 
Concurrent  JHKL  photometry  obtained  at  the  SAAO  was  used 
to  derive  Teg  based  on  the  integrated  bolometric  flux  and  the  fit¬ 
ted  Rosseland-mean  angular  diameter.  Parameters  of  the  best-fit 
model  atmospheres,  such  as  visual  phase,  effective  temperature, 
and  distances  are  consistent  with  independent  estimates,  which 
provides  additional  confidence  in  the  CODEX  modeling  approach. 

The  closure  phase  functions  of  our  targets  exhibit  non-zero 
values  at  all  wavelengths  with  a  wavelength  dependence  that  also 
correlates  with  the  positions  of  the  H2O  and  CO  bands.  This 


L7,  page  3  of  5 


A&A532,  L7  (2011) 


Table  3.  Results. 


R  Cnc 

XHya 

W  Vel 

RW  Vel 

‘/’Vis  (P) 

0.3  (356  d) 

0.7  (298  d) 

0.9  (397  d) 

(452  d) 

UD,  full  A  [mas]  (xl ) 

15.2  ±0.1  (17.2) 

7.8  ±0.1  (1.7) 

10.3  ±0.1  (4.0) 

14.0  ±0.1  (6.3) 

UD,  2.25  pm  (cont.)  [mas]  (xl ) 

13.2  ±0.1  (1.3) 

6.2  ±0.1  (0.3) 

8.0  ±  0.1  (0.2) 

10.6  ±0.1  (0.3) 

UD,  2.0  pm  (H20)  [mas]  (*;]) 

14.7  ±0.1  (4.8) 

8.3  ±0.1  (0.4) 

not  obtained 

not  obtained 

UD,  2.4  pm  (CO)  [mas]  (xl ) 

16.1  ±0.1  (8.6) 

9.0  ±0.1  (0.5) 

11.2  ±0.1  (1.0) 

14.6  ±0.1  (1.9) 

CP,  full  A  [deg] 

1 10  ±  4 

5  ±2 

-8  ±2 

14  ±3 

CP,  2.25  pm  (cont.)  [deg] 

45  ±5 

2  ±  2 

-4  ±  2 

3  ±2 

CP,  2.0  pm  (H20)  [deg] 

108  ±  4 

11  ±3 

not  obtained 

not  obtained 

CP,  2.4  pm  (CO)  [deg] 

153  ±4 

5  ±  3 

-12  ±  3 

34  ±4 

Best  P/M  Model  (xl ) 

Ml  In  (7.4) 

P40  (1.9) 

M25n  (4.8) 

M12  (9.5) 

Best  CODEX  Model  (xl) 

C50-378680  (4.6) 

O54-287940  (0.7) 

C8 1-250440  (1.6) 

O54-262360  (4.6) 

Oross  [mas] 

11.8  ±0.7 

4.9  ±  0.8 

5.8  ±  1.2 

9.2  ±  1.3 

rModel  [K] 

2628  ±  184 

3064  ±  380 

3216  ±  225 

3142  ±  243 

®Model 

0.45  ±  0.09 

0.67  ±  0.28 

0.77  ±0.12 

0.73  ±0.10 

Tegeas  [K] 

2694  ±  140 

2912  ±  310 

2990  ±  390 

mboi  not  avail. 

dpi  [pc] 

280  ±  28 

440  ±  44 

510  ±51 

366  ±  36 

4lodel  (©Ross,  -Csf)  [PC] 

252  ±  28 

465  ±  49 

362  ±  96 

217  ±  33 

t^Model  (Wtbol,  Tf^del)  [pc] 

246  ±  46 

503  ±  80 

400  ±  37 

mboi  not  avail. 

Notes.  The  visual  phase  ®ViS  and  period  P  was  obtained  from  fitting  a  sine  function  to  the  AAVSO  and  AFOEV  data  of  the  last  ten  periods.  For 
RW  Vel  such  data  are  not  available,  and  we  list  the  GCVS  (Samus  et  al.  2009)  period  instead.  The  UD  diameters  are  obtained  from  a  fit  to  all  data 
(full  A),  to  the  continuum  bandpass  2.24—2.26  /im  (2.25  pm),  to  the  H20  bandpass  1.95-2.05  pm  (2.0  pm),  and  to  the  CO  bandpass  2.35-2.45  pm 
(2.4 pm).  The  closure  phases  CP  were  averaged  over  the  same  bandpasses.  The  Rosseland  mean  angular  diameter  ©ross  is  an  average  of  the  fitted 
diameters  based  on  several  closely  fitting  CODEX  model  atmospheres  of  any  series  with  average  effective  temperature  r^odei  antj  an  vjsua|  phase 
®Modei-  T™™  was  derived  from  0ross  and  mboi  (Table  2).  d?L  is  the  period-luminosity  distance  from  Whitelock  et  al.  (2008).  RW  Vel  is  based  on 
the  same  scale,  but  using  the  GCVS  period  and  the  mean  K  magnitude  from  Smith  et  al.  (2002).  The  last  two  rows  give  the  distances  based  on 
@ross  and  Rr^1  of  the  models  at  the  respective  phase,  and  based  on  mboi  and  M^°del. 


result  indicates  a  complex  non-spherical  stratification  of  the  ex¬ 
tended  atmosphere  of  Mira  variables.  The  most  significant  devi¬ 
ation  from  point  symmetry  is  observed  in  the  H2O  band  around 
2.0 pm,  in  particular  for  the  clearly  resolved  target  R  Cnc.  We  in¬ 
terpret  this  signal  as  an  indication  of  inhomogeneities  or  clumps 
within  the  water  vapor  layer  on  scales  of  a  few  resolution  ele¬ 
ments  across  the  stellar  disk.  These  inhomogeneous  water  shells 
have  also  been  detected  for  the  Mira  variable  U  Ori  (Pluzhnik 
et  al.  2009)  and  the  symbiotic  Mira  variable  R  Aqr  (Ragland 
et  al.  2008).  The  deviations  from  point  symmetry  at  the  near¬ 
continuum  bandpass  may  be  related  either  to  photospheric  con¬ 
vection  cells  (cf.  Freytag  &  Hofner  2008)  or  the  inhomogeneities 
of  molecular  layers  that  may  also  contaminate  this  bandpass. 

Inhomogeneous  or  clumpy  molecular  layers  may  be  ex¬ 
pected  as  the  result  of  chaotic  motion  induced  by  the  interac¬ 
tion  of  pulsation  and  shock  fronts  with  the  extended  atmosphere. 
The  CODEX  models  indicate  that  the  outer  mass  zones  outward 
of  1 .5-2  Rosseland  radii  are  only  loosely  connected  to  the  stel¬ 
lar  pulsation  (cf.  Fig.  1  of  Ireland  et  al.  2008).  Likewise,  outer 
mass  zones  on  different  sides  of  the  star  may  be  only  weakly 
correlated  with  each  other  and  may  have  different  extensions. 
Icke  et  al.  (1992)  described  as  well  that  the  outer  layers  of  an 
evolved  AGB  star  may  respond  with  chaotic  motion  to  the  pul¬ 
sations  that  originate  in  the  stellar  interior.  Observations  of  wa¬ 
ter  vapor  layers  may  be  particularly  sensitive  to  these  effects,  but 
other  molecules,  such  as  CO,  that  are  expected  to  be  plentiful 
in  the  shocked  region  of  the  atmosphere  (cf.  CherchnefF  2006) 
would  also  be  affected  by  this  large-scale  chaotic  motion.  The 
resulting  clumpy  structure  of  molecular  layers  may  explain  our 
complex  wavelength-dependent  closure  phase  signal. 


Further  interferometric  campaigns  with  high  spatial  and 
spectral  resolution  are  clearly  needed  to  characterize  in  detail  the 
morphology  of  atmospheric  molecular  layers  in  Mira  variables. 
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Fig.  2.  Flux,  squared  visibility  amplitude,  closure  phase,  and  uniform  disk  diameter  (from  top  to  bottom )  as  a  function  of  wavelength  for  the 
Mira  variables  X  Hya,  W  Vel,  and  RW  Vel  (from  left  to  right).  As  Fig.  1,  but  for  the  remaining  targets  of  our  sample. 
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